Signature of ballistic effects in disordered conductors 
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Statistical properties of energy levels, wave functions and quantum-mechanical matrix elements 
in disordered conductors are usually calculated assuming diffusive electron dynamics. Mirlin has 
pointed out [Phys. Rep. 326, 259 (2000)] that ballistic effects may, under certain circumstances, 
dominate diffusive contributions. We study the influence of such ballistic effects on the statistical 
properties of wave functions in quasi-one dimensional disordered conductors. Our results support 
the view that ballistic effects can be significant in these systems. 



One of the simplest and yet most widely used mod- 
els for a disordered wire is that of independent 
electrons moving in a disordered potential described by 
a 5-correlated Gaussian random function 'L'(r) (where 
r = (x,y, z) is a spatial coordinate) 



(^(r)) =0, 
(u(r)w(r')) = (27ri/T)-i(5(r - r') . 



(1) 



Here (• • • ) is an average over potential realisations, v = 
1/{VA) the electronic density of states, V the volume of 
the system, A the mean level spacing, and t the scatter- 
ing time. The corresponding Hamiltonian is 



H 



P 

2m 



v{r) 
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with mass m and momentum p. For the model (jljL, 
it has been proven that on small energy scales (of 
the order of A) and for sufhciently weak disorder the 
quantum-mechanical correlation functions are universal 
and equivalent to those of random matrix theory (RMT) 
Moreover, deviations from universal statistics are 
parameterised by the dimensionless conductance g ~ 
2TThD/{L^A) where D is the diffusion constant and L 
is the linear extension of the system. In the limit of 
g — > oo, fluctuations are of RMT type. Deviations may 
be calculated perturbatively, as an expansion in g~^ (see 
for instance and references therein). 

These important results were established using an ef- 
fective field theory, the nonlinear a model (NLSM) 
valid for <^ £ (where Ap is the Fermi wave length) and 
on length scales larger than the mean free path £ = v-pT 
(where vp is the Fermi velocity). Efetov's NLSM as- 
sumes that the electrons undergo diffusive motion and, 
within a saddle-point approximation, statistical proper- 
ties of eigenvalues, wave functions, and matrix elements 
are characterised by a classical return probability ex- 
pressed in terms of the diffusion propagator. 

It has, however, been suggested that under certain cir- 
cumstances, ballistic trajectories (i.e., trajectories scat- 
tering just once before returning to their starting point) 



may influence the fluctuations in the model (|l|g) ||3| . At 
present it is unclear to which extent the fluctuations in 
(|l|,||) may differ from those obtained from the diffusive 
NLSM. It is thus necessary to quantify the significance 
of ballistic effects using alternative methods. 

One possibility would be to derive a more general 
NLSM for (0j^) incorporating ballistic effects ||] . Very re- 
cently, for instance, a ballistic NLSM was used to discuss 
statistical properties of energy levels and wave functions 
in ballistic quantum billiards with surface scattering Q. 

In this paper, we determine the statistical properties of 
wave-function amplitudes in the model (|l|,||) using exact- 
diagonalisation calculations for one- and quasi-one di- 
mensional (Id and quasi-lc?) tight-binding models with 
random on-site potentials. We determine to which ex- 
tent the statistics are influenced by ballistic effects. In 
addition to spatially uncorrelated disordered potentials 
[corresponding to (|l|)] we also consider potentials with 
smoothly varying spatial correlations (Eq. (|^) below) . In 
the latter case one expects ballistic contributions to be 
suppressed. Our results support the view that ballistic 
effects can be significant in quasi-one dimensional disor- 
dered systems. 

Formulation of the problem. Within the diffusive 
NLSM, one finds that deviations from universal statistics 
are parameterised by the time-integrated return proba- 
bility |,| 

in quasi-ld , (3a) 
\og{L/£) in 2d, (36) 
L/£ in 3d (3c) 



P{r,r-u^O)^ g- 



where P{r,r';uj) is the diffusion propagator. It was 
pointed out in |3| that the time-integrated return proba- 
bility may have additional contributions of the form 



P(r, 



0) 



g-^logi£/Xp) 
Af/^ 



in 2d, 
in 3d . 



(4a) 
(46) 



These arise from ballistic trajectories contributing to the 
return probability on small length scales. In 3d, the dif- 
fusive contribution behaves as ^ (Ap/^)^ which is much 
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smaller than (46). The applicability of the diffusive 
NLSM is thus questionable, an issue raised in [||J^ . In 2d 
systems, the effects are less drastic but can nevertheless 
be appreciable. In quasi-ld wires, the ballistic contri- 
bution is of the form (^) (when the sample is locally 
3d). Then the ratio of (3a) and (46) determines whether 
ballistic effects are important or not. 

The ballistic contributions (4a, 6) are suppressed in the 
case of a smoothly correlated disordered potential, such 



as 



(u(r)r;(r')) = (27ri.r)-i/(|r-r'|/4) 



(5) 



where f{x) is a dimensionless smooth function decaying 
on the scale a; ~ 1. When the correlation length is much 
larger than the Fermi wave length, ^ Ap, ballistic 
contributions of the type (jj) are negligible. 

The following questions arise: How important are bal- 
listic contributions of the type (^, 6) in low-dimensional 
disordered quantum systems? How exactly are they sup- 
pressed when £c is increased? Under which circumstances 
is the diffusive NLSM applicable? 

Method. In order to answer these questions we 
have performed exact diagonalisation studies of a tight- 
binding version of (|^) on a cubic lattice with lattice spac- 
ing ao 



H 



rx' 



r 



t 
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Here cj. and cr are the creation and annihilation opera- 
tors, the hopping amplitudes are |trr'| = 1 for nearest- 
neighbour sites and zero otherwise. We consider 6- 
correlated potentials 



[vrvr' 



rr' 



(7) 



as well as spatially smooth potentials 

{vrvr') = (W^Vl2)exph|r-rf /(24)']. (8) 

As usual, the parameter W characterises the strength 
of the disorder. As is well known, the eigenvalues Ej 
and wave functions ■0j(r) of this Hamiltonian, for d > 2, 
on small energy scales and for sufficiently weak disorder 
{g oo), exhibit fluctuations described by RMT. De- 
pending on the phases of the hopping amplitudes irr', 
Dyson's Gaussian orthogonal or unitary ensembles are 
appropriate Q. We refer to these cases by assigning, as 
usual, the parameter P — 1 to the former and /3 = 2 to 
the latter. 

By diagonalising the Hamiltonian (|^) using a modified 
Lanczos algorithm pO| , we have determined the distribu- 
tion function 

fpiE,r;t) = A{j2dit-\^P,{r)\^V)g,,iE-E,)) (9) 



of wave-function amplitudes for Id and quasi-ld metal- 
lic samples. The wave functions are normalised so that 
(IV'j(r)P) = V^^ and gr^{z) is a window function of width 
?7, centered around z = and normalised to unity. 

Results. We first discuss results for a chain of length 
L = 2 X lO'^ao with periodic boundary conditions, using 
the smoothly correlated potential (||). Fig. ^(a) shows the 
x-averaged distribution functions f{E;t) = {J{E,x]t))x 
as a function of t for M^ = 0.5, E = —1 and for several 
different values of £c. Our numerical results are fitted by 
the expression 



/(i.;t).lexp(-^ 



(10) 



derived in ||Tl[ for a Id chain of length L ^ ^. Here f 
is the localisation length. Eq. ( |lo|) describes the exact- 
diagonalisation results well for t ^ L/^ exp(— L/^), as 
expected |12 . Fig. |l|(b) shows ^ as a function of 4- We 
observe that ^ increases exponentially fast with £c for 
£c not too large. This is in keeping with second-order 
perturbation theory [O , valid for small W where 



oo 



96aal~E^/4: ^ 



-x^al/{2£^)^^2ixk 



(11) 



with E = 2cos(fc). Evaluating the sum over x using 
Poisson summation one obtains for E = ~1 and large 4 



36 ag 



V^£,exp[-4TT^£l/{9al)]. (12) 



We observe reasonable agreement between Eqs. (|Tl|) (and 
( |T^ ) for large £c) and the results of our simulations. For 
stronger disorder and shorter chains we observe signifi- 
cant deviations from (^l]) (not shown). 

We now discuss our results for quasi- Id wires. As is 
well known, in the limit of large dimensionless conduc- 
tance, the distribution function (^ tends to the RMT 
result 



_ / exp(- 
exp(- 



(13) 



The P = 1-distribution is often referred to as the Porter- 
Thomas distribution [|| . Within the diffusive NLSM one 
obtains for the s-averaged relative deviations 5fp{E] t) = 
{fp{E,x-t))Jff\t)-l 



5fp{E-t)^P 



3/4-3</2+t74 for/3^1, 
l-2t+t^/2 for/3 = 2 



(14) 



(see 1^ and references therein). Here P = {P{x,x;uj — 
0))x- Eq. (^ is expected to be valid for small values 
oi t {t < P^"^). We emphasise that, according to the 
diffusive NLSM, P is the same for /3 = 1 and 2: it is 
related to the diffusion propagator, a classical quantity. 
Interference effects are contained in the coefficients of the 
polynomials in t. 
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According to Refs. HQ, the diffusive NLSM is ex- 
pected to be valid under the following conditions 



1 < ^/Af < S/Xl 



(15) 



where S — Muq is the cross-section of the wire. 

How may ballistic contributions, such as those given by 
Eq. (Q) affect the prediction (p^)? The heuristic argu- 
ment of Ref. H] proceeds as follows: in the (5-correlated 
case, short ballistic trajectories contributing to the re- 
turn probability scatter only once before returning and 
are thus self-retracing. When fitting P independently for 
(3=1 and 2 to Eq. (|^, the values may differ. In the 
limit where ballistic effects dominate, P1/P2 — > 1/2. If 
ballistic effects are absent (for sufficiently smoothly cor- 
related Vy) one expects P1/P2 = 1- 

Comparing expressions (^) and (^), one sees that in 
the (5-correlated case the diffusive contribution is domi- 
nant only for small enough values of 5, i.e., when g <C 
£/Xf- Since g = ^/L ~ S£/{XpL), this implies a condi- 
tion 



or, if Ap — ao. 



S/Xl < X/Af 



M < L/aa . 



(16) 



(17) 



Therefore, we expect the ballistic effects to dominate for 
large M and small ic- 

We have evaluated 6fp{E,t) for (5-corrclated as well 
as for smooth disordered potentials in a L x 
sample of length L and cross section S. Our results for 
the (^-correlated case are summarised in Fig. |. The re- 
sults for the smooth case {£c/ao = 2) are similar and 
are shown in Fig. ||. In both cases, we have also plotted 
fits of Eq. (14) to the data, using the parameters Pi and 
P2 as a fitting parameters. We see from Figs. || and § 
that the functional forms of the deviations dfi3{E;t) are 
wefi described by (|l|). In Fig. ^ the values of P1/P2 
are shown, as obtained from fitting ( p| ) to the data in 
Figs. H and |[ for /3 = 1 and 2. Despite the fact that the 
scatter is large, we identify the following trends. 

First, in the smoothly correlated case, P1/P2 — 1 for 
all values of M. This corresponds to the result derived 
using the diffusive NLSM. In the case of a smooth poten- 
tial, large-angle scattering of the electrons is reduced and 
the contribution of ballistic, self-retracing trajectories to 
the time-integrated return probability is negligible. 

Second, in the (5-correlated case we observe a crossover 
from P1/P2 ~ 1 for smaU values of M to P1/P2 ~ 1/2 
for large values of M. The crossover takes place at values 
of M of the order of L/qq. This is consistent with the 
heuristic prediction ( pj| ) and indicates that ballistic con- 
tributions to the return probability are of the form (46) . 
They dominate the diffusive contribution (3a) for large 
values of M and are negligible for small values of M. In 
the latter case the diffusive NLSM applies, in the former 
case it does not. 



Conclusions. Using exact-diagonalisation calculations 
of the quasi-ld tight-binding model (6-8), we have deter- 
mined statistical properties of wave-function amplitudes 
and compared the results to predictions of the diffusive 
NLSM. Our results depend on the nature of the spatial 
correlations of the disordered potential which may be S- 
correlated or spatially smooth. The conclusions can be 
summarised in three points. First, in the (5-correlated 
case the diffusive NLSM applies provided the wires are 
not too thick (M ^ L/ao). This is consistent with the 
following heuristic ansatz: The classical return proba- 
bility is assumed to be a sum of diffusive and ballistic 
contributions If M ^ L/ao, the diffusive contri- 

bution (3a) dominates and the diffusive NLSM is ap- 
plicable. If M ^ L/ao, the ballistic contribution (46) 
dominates and the diffusive NLSM is not applicable (we 
find P1/P2 ~ 1/2 and not 1 as predicted by the diffusive 
NLSM). Second, the functional form of the deviations, is 
still very well approximated by the diffusive NLSM pre- 
diction even for thick wires. Third, in the smoothly cor- 
related case, the exact-diagonalisation results agree with 
the NLSM predictions. 

Our results confirm a heuristic discussion of ballistic 
effects 1^. Their infiuence on the distribution of wave- 
function amplitudes is thus qualitatively understood, but 
a quantitative understanding has not yet been achieved. 
It would be of great interest to calculate the distribution 
function fp{E,x;t) using the ballistic NLSM discussed 
in 0, and to compare to the results summarised in the 
present paper. 
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FIG. 2. Shows 5ff3{E;t) for the model (7), for = 1 (•) 
and /3 = 2 (□), for W ~ 0.7, L = 128 ao, and for different 
values of M. Also shown are fits to the data according to ( |l4| ) 
(lines), with P as a fitting parameter. 
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FIG. 1. (a) The distribution functions f{E;t) for a Id 
chain of length L/ao = 2 x 10*, with W = 0.5, for different 
values of Ic/ao (symbols). They have been fitted to Eq. ( [ic| ) 
(lines), using ^/ao as a fitting parameter, (b) The correspond- 
ing localisation lengths ^ as a function of i^. Also shown are 
expressions (0) (full line) and (^ (dashed line). 
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FIG. 3. The same as Fig. | but for the model (8), for 
/3 = 1 (•) and /3 = 2 (□), 4/ao = 2.0, and W = 0.5. 
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FIG. 4. P1/P2 versus VM for 5-correlated (□) and 
smooth (o) potentials. The parameters are as in Figs. ^ 
and M. 
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